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ABSTRACT 


^  The  pick-up  and  delivery  problem  (PDP)  is  a  special  case  of  the 
vehicle  routing  problem  and  can  be  defined  as  follows.  There  exists  N 
tasks  to  be  performed  by  a  fleet  of  vehicles.  Each  task  consists  of 
picking-up  a  certain  amount  of  goods  from  its  supply  node  and 
transporting  these  to  its  associated  demand  node.  The  problem  is  to 
perform  these  N  tasks  so  that  a  certain  objective  is  optimized.  The  key 
characteristics  of  the  PDP  are  the  precedence  relations  in  visiting  the 
supply-demand  nodes  and  fluctuations  in  vehicle  capacity  usage  during  the 
course  of  the  tour. 


Two  linear  binary  models  which  incorporate  the  special  PDP 
characteristics  are  presented.  The  first  model  is  for  the  single  vehicle 
problem.  This  model  is  extended  to  the  multiple  vehicle  case.  The 
relaxation  of  certain  constraints  in  the  single  vehicle  PDP  model  yields 
lower  bounds  for  the  optimal  solution.  These  lower  bounding  problems  are 
solved  by  dynamic  programming. 


The  lower  bounds  obtained  are  implemented  in  a  branch  and  bound 
procedure.  The  algorithm  is  able  to  solve  an  eight  task,  single  vehicle 
problem  optimally.  By  dualizing  one  of  the  PDP  constraints  and  solving 
the  resulting  program  by  a  Lagrangian  ascent  procedure,  better  quality 
lower  bounds  are  obtained.  This  permits  twelve  task,  single  vehicle 
problems  to  be  solved  optimally. 

Key  words:  Vehicle  Routing,  Pick-up  and  Delivery  Problem,  Dynamic 
Programing,  Lagrangian  Relaxation,  Branch  and  Bound 
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1.  INTRODUCTION 


Vehicle  routing  and  scheduling  problems  have  attracted  considerable 
attention  among  researchers.  Vehicle  routing  and  scheduling  is  a  term 
given  to  a  wide  class  of  problems.  Each  problem,  involves  routing  and 
scheduling  of  vehicles  operating  in  different  environments  and  performing 
different  tasks.  Bodin  and  Golden  [1]  attempt  to  classify  vehicle 
routing  and  scheduling  problems  into  different  groups.  Each  group 
concentrates  on  one  aspect  of  a  particular  problem.  One  group  of 
problems  is  identified  as  "pick-up  and  delivery  problems,"  (PDP).  In 
this  paper  pick-up  and  delivery  problems  are  analyzed  and  solution 
strategies  and  procedures  for  the  single  vehicle  case  are  proposed. 

1.1  Pick-up  and  Delivery  Problems 

PDP  can  be  described  as  follows: 

There  exist  N  tasks  to  be  performed  with  a  fleet  of  V  vehicles. 
Each  task,  r  -  1,...,N,  consists  of  picking-up  qf  units  from  a 
source  node  a(r)  and  transporting  them  to  the  corresponding  demand 
node  b(r).  Each  task  can  be  specified  by  three  parameters: 

(a(r) ,b(r) ,qf) .  The  problem  is  to  perform  these  N  tasks  with  the 
available  vehicle  fleet  so  that  a  certain  objective  is  optimized. 

PDP  can  be  classified  with  respect  to  different  schemes.  Bodin  and 
Golden  suggest  several  classifications  including: 

A.  Fleet  Size: 

1.  Single  vehicle 

2.  Multiple  vehicle 


B.  Operation  Type: 


1.  One- to-Many-to-One 

2.  Many-to-Many 

3.  Any  combination  of  these 

C.  Type  of  Fleet: 

1.  Homogeneous  case  (all  vehicles  are  identical) 

2.  Heterogeneous  case  (not  all  vehicles  are  identical) 

D.  Objective: 

1.  Minimize  routing  costs 

2.  Minimize  sua  of  fixed  and  variable  costs 

3.  Minimize  fleet  size 

E.  Time  Constraints: 

1 .  On  routes 

2.  On  tasks 

3 .  Mixed 

F.  Other  (problem  dependent)  constraints: 

1.  On  vehicle  loading/ unloading 

2.  On  tasks 

3.  On  operations 

The  nature  of  tasks  and  operation  types  are  the  two  main 
characteristics  of  PDPs  which  differentiate  them  from  the  general  vehicle 
routing  problems  (VRPs)  described  by  Christofides  [14].  In  PDPs,  both 
pick-ups  and  deliveries  are  Involved  where  in  VRPs  only  pick-ups  or  only 
deliveries  are  performed.  Because  of  this  added  complexity  of  PDPs,  the 
solution  strategies  and  procedures  which  are  effective  for  VRPs  cannot  be 
used  directly.  The  main  reasons  are: 

(1)  In  VRPs,  once  tasks  are  assigned  to  vehicles,  the  route  for 
each  vehicle  can  be  determined  by  the  solution  of  a  Travelling  Salesman 
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Problem  (TSP).  In  PDPs,  even  If  Casks  can  be  aslgned  Co  Che  vehicles, 

Che  routing  of  a  parclcular  vehicle  is  a  harder  problem  chan  a  TSP,  due 
Co  Che  precedence  reladons  of  Che  Casks . 

(2)  In  VRPs,  Che  capacity  usage  of  Che  vehicle  is  deCermined  aC  Che 
incial  node  of  ics  rouCe  and  Che  capaciCy  usage  declines  sCeadily  (in  Che 
delivery  case)  ChroughouC  Che  rouCe.  In  Che  case  of  PDPs,  ic  is  noc 
possible  Co  deCermine  Che  capaciCy  usage  in  advance  since  ic  fluctuates 
as  deliveries  and  pick-ups  are  performed. 

(3)  Operation  type  is  anocher  factor  which  differentiate  VRPs  from 
Che  PDPs.  In  VRPs,  operations  can  be  identified  as  one-Co^many  for 
single  domicile  problems  and  fev-to-many  for  the  multiple  domicile 
problems,  where  a  domicile  is  Che  node  at  which  goods  are  initially 
scored  (in  delivery  case).  In  both  of  these  problems  Che  goods  Co  be 
transported  can  be  considered  as  homogeneous  in  the  sense  that  they  are 
supplied  from  specific  domiciles.  In  the  case  of  PDPs  Che  goods  can  not 
be  considereed  as  homogeneous,  because  each  demand  point  of  a  Cask  can 
only  be  supplied  from  its  respective  supply  point.  This  non-homogenity 
of  the  goods  creates  a  multi-product  environment  which  further 
complicates  Che  problem. 

In  this  paper  "many-to-many"  PDPs  will  be  considered. 

"many-to-many"  means  that  the  pick-up  and  delivery  points  of  each  task 
are  all  distinct  points.  In  one( few)-to-many-to-one( few)  PDPs,  vehicles 
pick-up  goods  from  the  domiciles  and  delivery  them  to  "many”  demand 
points  and  then  pick-up  some  goods  from  these  demand  points  and  deliver 
them  to  the  domiciles.  Aa  far  as  the  operations  are  concerned 
"many-to-many"  problems  can  be  considered  as  a  generalization  of  the 
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"one-to-many-to-one"  problems.  The  solution  procedures  for  many-to-many 
problems  can  easily  be  adpated  to  the  later  problem. 

Pick-up  and  delivery  problems  are  among  the  most  frequently  faced 
problems  In  transportation  and  distribution.  In  practice,  these  problems 
arise  In  different  transportation  fields,  such  as  naval,  air,  rail,  and 
land  transportation.  Also,  these  problems  occur  In  different  environ¬ 
ments  such  as  urban  transportation.  Interstate,  and  Intercontinental 
transportation.  Some  reported  applications  of  the  PDPs  are: 

a)  Tractor  trailer  routing  and  scheduling;  [5],  [7],  [11] 

b)  Navy  Supply  Center  routing:  [22],  [29] 

c)  Static  dial- a- ride:  [10],  [12],  [27] 

d)  School  bus  routing  and  scheduling:  [2],  [20],  [24],  [28] 

e)  Air  Force  logistics:  [1] 

f)  Tanker-vessel  scheduling:  [3],  [9],  [19]. 

1 .2  Survey  of  Solution  Methods  for  the  Pick-up  and  Delivery  Problems 

Unlike  VRPs,  there  exists  little  or  no  general  literature  concen¬ 
trating  on  PDPs.  Each  of  the  studies  involves  a  particular  PDP  and  pro¬ 
poses  solution  procedures  which  makes  use  of  the  special  characteristics 
of  the  particular  problem  addressed.  Also,  there  currently  exists  no 
exact  procedures(s)  which  solves  these  problems.  According  to  Psaraftls 
[25],  even  moderate  size  static  dlal-a-ride  problems,  (N  <  10)  cannot  be 
solved  exactly  with  a  reasonable  computational  effort.  Most  of  the 
solution  methods  are  Inexact  procedures. 

1.2.1  Full  I/} ad  PDP 

The  PDPs,  where  •  C  ■  1  for  all  tasks,  are  called  the  full  load 
PDPs.  In  the  absence  of  time  restrictions,  the  full  load  pick-up  and 
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delivery  problem  is  formulated  by  Bodin  et  al  in  two  different  equally 
effective  ways.  One  approach  is  based  on  the  travelling  salesman 
formulation,  where  the  vehicle  is  restricted  to  visit  the  delivery  node 
of  a  particular  task  Immediately  after  visiting  its  associated  pick-up 
node.  The  same  problem  can  also  be  formulated  as  a  Chinese  postman 
problem  (CPP),  where  the  arcs  to  be  visited  are  identified  as  the  arcs 
between  the  pick-up  and  delivery  points  of  each  task.  Depending  on  the 
formulation  any  TSP  or  CCP  procedure  can  be  applied  to  the  problem. 

For  time  constraints.  Ball  et  al  [7]  propose  two  different 
heuristics.  The  first  is  a  two  phase  method  where  the  first  phase  con¬ 
sists  of  using  either  a  TSP  or  a  CPP  procedure  to  create  a  large  route. 

In  the  second  phase,  this  large  route  is  partitioned  into  smaller  feasi¬ 
ble  routes  which  can  be  handled  by  the  available  vehicles.  In  the  second 
method,  a  greedy  insertion  routine  is  used  in  which  origin/destination 
pairs  are  added  to  the  vehicle  routes  iteratively. 

1.2.2  Partial  Load  PDP 

PDPs  which  have  task  loads  less  than  the  vehicle  capacities  are 
called  partial  load  problems.  The  partial  load  problem  can  also  solve 
full  load  or  static  dial-a-ride  problems.  School  bus  routing,  airplane 
routing,  and  ship  scheduling  all  fall  under  partial  load  problems.  All 
procedures  available  in  the  literature  attempt  to  solve  this  problem 
heuristically . 

All  et  al  [1]  generate  "potentially  good  feasible"  routes  and,  with 
a  mixed  integer  program,  try  to  select  a  set  of  best  routes  satisfying 
routing  and  scheduling  requirements.  The  authors  apply  the  procedure  to 
a  60  base,  61  task  problem  of  the  U.S.  Air  Force.  As  they  point  out,  the 


selection  of  "potentially  good"  routes  is  the  crucial  step.  As  the 
ability  to  identify  "good"  routes  increases,  the  quality  of  the  resulting 
solution  also  Improves. 

The  second  approach  for  these  type  of  problems  are  the  TSP  based 
heuristics.  Inexact  TSP  procedures  are  modified  with  respect  to  the 
special  characteristics  of  each  problem. 

Assad  et  al  [5]  propose  a  two  phase  approach  to  a  partial  load  PDP. 
Phase  one  contains  the  routing  component  and  the  second  phase  involves 
the  scheduling  of  the  routes  found  in  the  first  phase.  The  routing  phase 
involves  an  insertion  algorithm  where  origin  and  destination  nodes  of  a 
task  are  inserted  into  one  of  the  existing  vehicle  routes.  At  the  end 
of  the  routing  phase  the  vehicle  routes  are  feasible  with  respect  to: 

1)  precedence  relations  for  each  task, 

2)  vehicle  capacities,  and 

3)  maximum  route  lengths. 

The  scheduling  phase  scans  each  vehicle  route  from  start  to  finish 
and  attaches  times  to  each  of  the  operations  involved.  This  phase 
minimizes  the  service  times  of  each  vehicle.  If  any  infeasibilities  are 
found  with  respect  to  time  constraints  they  are  ignored. 

Whinchel  et  al  [29]  propose  a  single  phase  method  for  the  scheduling 
of  palletized  cargo  delivery  among  the  warehouses  in  a  Navy  Supply 
Center.  In  their  approach,  vehicle  routes  are  created,  one-at-a-time ,  by 
attaching  a  task  to  the  end  of  the  route  under  construction.  The 
criterion  used  is  the  savings  concept  of  Clark  and  Wright  [18].  For  the 
same  problem,  an  insertion  routine,  which  examines  all  possible 
insertions,  was  developed  by  Jarvis  and  Kirca  [22].  This  procedure 
resulted  in  improved  solutions  of  the  problem. 
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1.2.3  Static  Dial-a-Ride 

The  static  dial-a-ride  problem  is  one  of  the  most  widely  studied 
pick-up  and  delivery  problems.  Procedures  for  the  dial-a-ride  problem 
will  be  examined  under  two  parts:  single  vehicle  and  multiple  vehicle 
cases . 

1.2. 3.1  Single  Vehicle  Dial-a-Ride  Problem 

When  there  are  no  time  windows  on  the  tasks,  Psaraftls  [25]  proposes 
an  optimal  procedure.  He  presents  a  mathematical  program  for  this 
problem  and  solves  it  by  dynamic  programming.  The  major  disadvantage  of 
his  approach  is  the  "curse  of  dimensionality''.  There  exist  N  state 
variables  in  each  of  the  2N  stages  in  his  formulation.  The  required 
computer  memory  is  on  the  order  of  3N.  Due  to  this  high  memory 
requirements  his  approach  fails  for  most  applications.  Computational 
experience  with  his  method  did  not  extend  beyond  problems  having  nine 
tasks.  Even  eight  and  nine  task  (customer)  problems  required  149  and  591 
CPU  seconds  respectively  in  finding  the  optimal  solutons. 

Psaraftls  [26]  also  develops  two  heuristic  proedures  for  the  same 
problem.  In  the  firt  heuristic,  an  approximate  TSP  tour  is  constructed 
using  Christofides'  [13]  procedure.  Then,  employing  improvement 
routines,  the  initial  TSP  tour  is  forced  to  be  feasible  with  respect  to 
both  precedence  and  capacity  constraints.  The  second  heuristic  uses  the 
minimal  spanning  tree  portion  of  Christofides'  algorithm.  The  initial 
TSP  tour  is  formed  and  then  this  tour  is  improved.  Psaraftls  also  shows 
that  his  first  algorithm  has  a  worst  case  bound  of  3.0. 

Armstrong  [4]  proposes  a  dynamic  programming  based  exact  algorithm 
for  the  problem.  Although  in  his  work  the  general  PDP  is  discussed,  all 
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the  procedures  developed  and  tested  are  for  the  static  dial-a-ride 
problem.  The  exact  algorithm  was  able  to  solve  highly  restricted  15 
customer  (task)  dial-a-ride  problems.  Armstrong  also  proposed  several 
heuristic  procedures  which  utilize  the  Clarke-Wright  savings  technique, 
greedy  approach,  and  an  insertion  routine,  he  attempted  to  improve  the 
tours  found  by  each  of  the  above  procedures  by  a  3-opt  improvement 
routine.  Computational  results  demonstrate  that  the  procedures  which 
selects  the  best  of  five  3-optimal  random  tour»  perform  better  than  the 
other  procedures  considered  in  his  work. 

For  the  case  involving  time  windows,  Sexton  and  Bodin  [27]  and  Baker 
[6]  propose  two  procedures.  Sexton  and  Bodin  [35]  try  to  minimize  the 
total  Inconvenience  which  a  customer  may  experience  due  to  excess  ride 
time  and  deviation  from  the  desired  delivery  time.  They  formulate  a 
mixed  integer  program  and  apply  Benders'  partitioning  method  to  this 
formulation.  However  the  partitioned  program  is  not  solved  exactly. 
Instead  they  make  use  of  the  master-subprogram  mechanism  of  this 
partitioning  in  order  to  find  a  heuristic  solution.  The  master  program 
is  a  routing  scheme  where  the  order  of  pick-ups  and  deliveries  are 
decided.  In  the  scheduling  phase,  the  exact  times  of  pick-ups  and 
deliveries  are  specified.  Then,  with  respect  to  the  schedules,  the 
objective  function  coefficients  of  the  routing  problem  are  adjusted. 

This  procedure  alternates  between  the  routing  and  scheduling  problems 
until  no  improving  solutions  can  be  found. 

Baker  [6]  proposes  another  method  which  minimizes  the  total 
servicing  time  of  the  vehicle,  he  formlates  a  nonlinear  program  and 
atempts  to  solve  this  program  by  several  procedures. 
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1.2. 3. 2  Multiple  Vehicle  Dial-a-Ride  Problem 

There  do  not  exist  effective  exact  solution  procedures  for  these 
type  of  problems.  Cullen,  Jarvis  and  Ratliff  [18]  and  Bodin  and  Sexton 
[10]  present  two  heuristic  procedures. 

Cullen  et  al  propose  an  interactive  optimization  scheme  to  develop 
“good”  vehicle  routes.  The  promising  alternatives  are  suggestd  through 
the  use  of  shadow  prices  for  a  set  partitioning  formulation  of  the 
problem. 

Bodin  and  Sexton  propose  a  three  step  approach.  In  the  first  step 
customers  are  assigned  to  vehicles.  The  second  step  attempts  to  improve 
this  assignment .  In  the  final  step  each  vehicle  route  is  developed .  As 
the  authors  point  out,  finding  an  initial  set  of  "good"  clusters  is  still 
an  unsolved  problem.  They  propose  two  alternative  methods  for  the 
initial  assignment  of  tasks  to  the  vehicles.  The  authors  applied  the 
algorithm  to  the  problem  of  a  subscriber  dial-a-  ride  system  and  were 
able  to  find  improved  solutions  compared  to  those  of  manual  scheduling. 


2.  MODELS  FOR  THE  PICK-UP  AND  DELIVERY  PROBLEM 


Christof ides  discusses  mathematical  modeling  approaches  for  vehicle 
routing  problems.  Although  his  discussion  is  for  the  general  VRP,  most 
of  his  comments  about  the  modeling  approaches  are  valid  for  PDPs.  He 
Identifies  two  schemes  for  the  routing  problem,  namely,  set  partitioning 
based  models  and  TSP  based  models.  Further,  he  discusses  some  advantages 
and  disadvantages  of  each  scheme. 

One  disadvantage  of  a  set  partitioning  model  is  that  for  most 
practical  problems  it  is  impossible  to  generate  and  store  all  the 
alternative  feasible  routes.  On  the  other  hand,  a  TSP  based  model  is  not 
able  to  express  most  realistic  constraints  explicitly. 

2.1  Single  Vehicle  Pick-up  and  Delivery  Model 

Consider  a  pick-up  and  delivery  problem  where  there  exist  N  tasks 
and  a  vehicle  with  a  capacity  of  C  units.  The  following  characteristics 
will  be  assumed  for  the  operation  of  the  system: 

1)  The  vehicle  is  initially  stationed  at  the  base  located  at 
node  1. 

2)  The  route  starts  and  ends  at  the  base. 

3)  No  task  splitting  is  allowed.  That  is,  if  a  supply  node  is 
visited  the  vehicle  picks-up  the  entire  load  and  delivers  it  to 
the  associated  demand  node  on  its  first  visit. 

4)  All  supply  and  demand  nodes  are  distinct  in  the  sense  that  each 
supply  node  supplies  only  one  demand  node  and  each  demand  node 
obtains  it's  requirements  from  only  one  supply  node.  (This 
results  in  no  loss  of  generality.) 


S)  The  objective  Is  to  service  all  tasks  with  a  minimum  route 
length. 

Let 

N  “  number  of  tasks  to  be  serviced 

M  -  total  number  of  nodes.  (M  -  2N  if  the  base  node  Is  also  a 
supply  node,  or  M  ■  2N  +  1  otherwise) 

S  ■  set  of  supply  nodes 

D  «  set  of  demand  nodes 

a(r)  -  supply  node  of  task  r 

b(r)  -  demand  node  of  task  r 

q  -  units  of  load  for  task  r 

C  -  capacity  of  the  vehicle 

d^j  *  distance  between  node  1  and  node  j 

In  the  classical  TSP  models,  based  on  the  assignment  of  arcs,  the 

decision  variable  determines  which  arcs  are  in  the  tour  and  which  arcs 
are  not  visited  by  the  tour.  In  such  a  formulation  the  sequence  of  node 
visits  are  not  readily  available,  and  the  model  does  not  keep  track  of 
the  sequence  of  nodes  visited.  In  the  case  of  PDPs,  In  order  to  enforce 
the  precedence  and  capacity  constraints,  it  is  necessary  to  keep  track  of 
the  sequence  of  node  visits  as  well  as  the  arcs  traversed.  The  model 
must  track  the  sequence  of  the  nodes  visited  so  that  a  visit  to  the 
demand  node  of  a  particular  task  can  be  restricted  to  a  position  higher 
in  the  sequence  than  its  associated  supply  node.  Also,  the  model  should 
be  able  to  keep  track  of  the  load  on  the  vehicle  so  that  the  vehicle 
capacity  constraint  can  be  enforced. 

If  there  are  N  tasks  to  be  performed,  the  vehicle  has  to  make  M  -  2N 
+  1  stops,  i.e.  one  stop  for  each  pick-up  or  delivery  node  plus  one  for 
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Che  base  node.  The  problem  is  to  determine  Che  nodes  visited  at  each 
stop,  so  that  all  tasks  are  performed  with  a  minimum  tour  length.  The 
problem  could  be  formulated  as  an  M-stage(stop)  problem,  where  the  node 
visited  at  each  stage  is  identified. 

The  model  presented  here  considers  the  routing  problem  as  a 
stage-by-stage  process.  In  each  stage  the  vehicle  moves  from  one  node  to 
another.  Since  there  are  M  nodes  to  be  visited,  the  model  makes  M 
decisions,  one  in  each  stage.  Let 


if  the  vehicle  moves  from  node  i  to  j  at  stage  t 
otherwise 


Then  the  single  vehicle  pick-up  and  delivery  model  becomes 


M  M  M 

Min  Z  Z  Z  d,  .x^, 
t-1  i-1  j-1  J  3 


.t. : 
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j-2 
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j-1 
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M  M 
Z  Z 
t-1  j-1 
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(2.1) 


(2.2) 


(2.3) 


t  -  2, . . . ,M  (2.4) 


i  -  2, . . . ,M  (2.5) 
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Z 

r-1 


M  t  M  t-l 

E  Xij  '  E  Xji  "  ° 

j-1  j-1 


t  M 

Z  Z  q„(x 
£-1 


A  VVrtJ  -  Vr)J>  ‘ 


M  M 


A  A  w >x 


■  0  or  1 


2 . .  i 


t 


r 


i.J.t 


(2.6) 

2 . M  (2.7) 

1,...,N  (2.8) 

(2.9) 


Constraints  (2.2)  and  (2.3)  state  that  the  route  must  start  and  end 

at  the  base  node  1.  Constraints  (2.4)  force  the  vehicle  to  visit  exactly 

one  node  at  each  stage  and  (2.3)  restricts  it  to  visit  each  node  exactly 

once.  Constraints  (2.6)  are  the  flow  equations  for  each  node  at  each 

stage,  staging  that  if  the  vehicle  enters  a  node  at  stage  (t-l)  then  it 

must  leave  that  particular  node  at  stage  t.  Constraints  (2.2)  through 

(2.6),  together  with  the  integrality,  constraint  (2.9),  define  a  TSP  tour 

which  begins  and  ends  at  node  1. 

Constraints  (2.7)  are  the  vehicle  capacity  constraints  for  each 

stage  and  constraints  (2.8)  are  the  task  precedence  constraints.  These 

restrict  the  vehicle  to  visit  a  supply  node  before  visiting  its 

3  2 

associated  demand  node.  The  model  has  M  binary  variables  and  M  +2+4M 
constraints. 

Time  window  constraints  on  the  tasks  can  be  imposed  on  the  model  as 
follows.  Let 


(e^,e^)  -  earliest  and  latest  pick-up  or  delivery  times  at  node  i 
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z 

t 


elapsed  time  at  stage  t.  ( z.  ■  0) 

time  to  move  from  node  1  to  j. 

total  loading  and  unloading  time  at  node  1 


Then, 


M  M 

z,.  .  +  E  1  (p,.  +  Ox*  -  z  t  -  2, . . . ,M  (2.10) 

C  j-1  1-1  J  J 


M  M  M  _  M 

E  £i  £  xA.  <  z  <  £  eA  £  x^.  t  -  2 . M  (2.11) 

1-1  1  j-1  C  1-1  1  j-1  J 


Note  that  (2.10)  defines  the  elapsed  time  at  stage  t  and  (2.11) 
enforces  the  time  windows  for  each  stage.  If  node  s  Is  visited  at  stage 
M  t  M  t 

t,  then  £  x  -  1  and  £  x,  -  0  for  all  1  -  s,  hence  (2.11)  becomes: 
j-l  8J  j-1 


<  z - 


<  e 


s 


Constraints  (2.10)  and  (2.11)  can  be  combined  into  a  single  equation 
of  the  form: 


for  t  -  2, . .  .,M 


M  M  t-1  M 

£  e  £  xj  <  £  £ 

1-1  j-1  J  W  1-1 


M 

£ 


j-1 


+  v^x 


l 

ij 


M 

<  £ 
1-1 


M 

£ 


j-1 


(2.12) 
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2.2  Multiple  Vehicle  Pick-up  and  Delivery  Model 

The  multiple  vehicle  case  can  be  formulated  as  an  extension  of  the 
single  vehicle  model.  Using  the  assumptions  of  Section  2.1,  suppose 
there  are  V  vehicles  available  to  service  N  tasks.  Let  the  capacity  of 
vehicle  k  be  C,^  respectively.  Further  let 


yrk 


1  if  task  r  is  processed  by  vehicle  k 

i 

0  otherwise 

■ 

■» 

1  if  vehicle  k  moves  from  node  i  to  j  at  stage  t 

< 

0  otherwise 


Then  the  model  becomes 

M  M  M  M 
Min  EE  IE  d, 

k-l  t-1  1-1  J-l  «  1Jk 

s.t. 
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t-1 
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t-1 
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j-l 
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Z 

j-l 


V 

£  y 

k-l 


rk 


a(r)jk 


*b( r) jk 


(2.13) 


(2.14) 


r  *  1c  *  2,  •  •  • , V 


(2.15) 


r  *  1 j • • • jNj  k  m  1 | | V 


M  M 

**  2  xiik  <  1  t  m  2****»M»  k  “  l»***»v  (2.16) 

i-1  j-l  J 


^  Xijk  "  Xjik  “  Q  1  "  2»  *  * •  *M*  c  "  2 . .  k  "  1 . V  (2.17) 
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M  1 

A " 1 


k  -  2, . . .  ,V  (2.18) 


t  N  H 
SEE 
1-2 


A  j-i  qr(Xa(r)jk”  Aojk5  <  Ck 


1  —  2| * • « |M 


k  - 


(2.19) 


M  M 


£1  j\  t<Xb<r)Jk  ~  VOjk5  >  1  r  *  k  -  (2.20) 


Xijk  "  0  °r  1  1»J»t  “  1 . M;  k  “  1 . V 


(2.21) 


yrk  “  0  or  1  r  “  !»•••.»;  k  *  1 . v  (2.22) 


In  constraints  (2.14)  each  task  is  assigned  to  one  vehicle. 
Constraints  (2.15)  ensure  that  the  supply  and  demand  nodes  of  each  task 
are  visited  only  by  the  vehicle  to  which  the  task  Is  assigned.  The 
remaining  constraints  ensure  feasible  routes  for  the  vehicles  similar  to 
the  single  vehicle  model.  This  model  has  V  times  as  many  binary 
variables  and  constraints  as  the  single  vehicle  model. 

Both  the  single  and  multiple  vehicle  models  are  large  binary 
problems  which  are  not  likely  to  be  solved  by  direct  optimization 
methods.  In  the  next  section  some  alternative  procedures,  based  on  the 
above  formulations,  will  be  discussed.  These  may  be  helpful  in  solving 
some  moderate  size  problems. 
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III.  LOWER  BOUNDS  FOR  THE  SINGLE  VEHICLE  PICK-UP  AND  DELIVERY  PROBLEM 
Christofldes  et  al  [15,16]  propose  a  class  of  lower  bounding 
procedures  for  the  TSP  and  VRPs .  They  are  based  on  state  space 
relaxations  of  the  dynamic  programming  formulation  of  the  respective 
problems.  Using  these  Ideas,  potentially  effective  lower  bounding 
schemes  can  be  developed  for  the  formulation  given  in  Section  2.1. 


3.1  Lower  Bounds  Through  Minimum  M- tours 

Consider  the  single  vehicle  model  of  Section  2.1.  Suppose  the 
constraints  (2.5),  (2.7),  and  (2.8)  are  relaxed.  Then  the  resulting 


program  is 


M  M  M 

Min  Z  Z  *■  d.  ,x* 
t-1  i»l  j-1 


(3.1) 


s  .t .  Z  x  ■  1 
j-2 


(3.2) 


(3.3) 


M  M 
I  E  x' 
i-1  j-1 


t  -  2, . . . ,M  (3.4) 


M  M 

£  x'  -  E  -  0 

j-1  J  j-1  3 


t  -  2, . . . ,M;  i  -  2, . . . ,M  (3.5) 


x!j  -  0  or  1 


i,i,t  -  1 , . . . ,M  (3.6) 
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In  the  above  program,  (3.1)  through  (3.6)  define  a  minimum  M-tour 
which  visits  M  nodes  and  begins  and  ends  at  node  1.  Since  constraints 


(2.5)  are  relaxed,  the  tour  may  contain  subtours  in  the  sense  that  a  node 
may  be  visited  more  than  once.  Further,  the  tour  may  not  be  feasible 
with  respect  to  the  capacity  and  precedence  constraints  of  (2.7)  and 
(2.8)  respectively.  As  Houck  et  al  [21]  show,  the  value  of  the  minimum 
M-tour  is  a  lower  bound  for  the  TSP  and  also  for  the  single  vehicle  PDP. 

The  minimum  M-tour  problem  is  a  shortest  path  problem  which  can  be 
solved  efficiently  by  dynamic  programming.  In  solving  the  problem,  2- 
node  sub tours,  which  are  of  the  form  i-j-i,  can  be  prohibited  by  the 
procedure  suggested  by  Christofides  [14]. 


3.2  Lower  Bounds  Through  Minimum  QM-to urs 


Let  be  the  total  number  of  units  of  load  that  should  be  serviced 
by  the  vehicle,  i.e. 


(3.7) 


We  can  add  the  following  additional  constraint  to  the  M-tour  model 
of  Section  3.1.1  (which  was  redundant  for  the  single  vehicle  model). 


M  N  M 

Z  1  Z  q  xC,  . 

t-l  r-l  J-l  r  *<r)J 


(3.8) 


Constraint  (3.8)  forces  the  minimum  M-tour  to  visit  a  combination  of 
supply  nodes  such  that  a  total  of  units  of  load  are  picked-up.  Such 
an  M-tour  will  be  called  a  QM-tour  and  the  value  of  this  tour  Is  a  valid 
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lower  bound  on  the  single  vehicle  model.  The  minimum  QM-tour  problem  can 
be  solved  by  dynamic  programming  as  follows.  Let 


f^CQ)  “  length  of  the  shortest  path  from  node  1  to  i  which  visits  t  nodes 
and  picks  up  a  total  load  of  Q. 


Then  for  i  “  2,...,M;  t  ■  3,...,M;  and  Q  ■  0,...,Q 


f.(Q)  “  Minimum 


(Q-q<)  +  d,.  if  its 


T  Qji 


fr1(Q>  +  ^ 


if  ieD 


where , 


fJ(Q)  - 


d^  if  Q  -  q1  and  i  e  S 


otherwise 


Then  the  length  of  the  minimum  QM-tour  can  be  computed  as  follows. 


LB 2  “  Minimum  {fJ(Qx)  +  ) 


3.3  Lower  Bounds  Through  Minimum  CM-tours 

The  capacity  constraints  (2.19)  of  the  single  vehicle  model  can  be 
Imposed  on  the  M-tours  such  that  the  vehicle  capacity  would  not  be 
violated  at  any  stage.  The  minimum  M-tour  which  does  not  violate  the 
capacity  constraints  will  be  called  a  CM-tour  and  the  dynamic  programming 
recursive  equations  can  be  written  as  follows.  Let 
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f  ^(U)  “  length  of  the  shortest  path  from  node  1  to  i  visiting  t  nodes 
and  having  a  total  load  of  U  units  at  node  i. 


Then  for  all  i  *  t  ■  3,...,M  and  U  »  0,...,C 


f  .  (U)  ■  Minimum 


fr1<o-v + dn 

fr1(° + "i> + dji 


i  e  S 


i  e  D 


where 


fj(U) 


d,,  if  U  -  q  and  i  e  S 
li  ^r 


otherwise 


Then  the  length  of  the  minimun  CM-tour  is  a  lower  bound  on  the 
single  vehicle  model ,  and  is  given  by 


LB3  ■  Minimum  (f^(0)  +  } 

ieD 


At  stage  M,  the  load  at  every  node  must  be  zero  since  all  the 
deliveries  must  be  completed  up  to  this  stage. 

3.4  Lower  Bounds  Through  Minimum  CQM- tours 

Both  the  capacity  and  the  total  load  constraints  (3.8)  can  be 
enforced  on  the  minimum  M-tours.  Such  an  M-tour  will  be  called  a  CQM- 
tour  and  the  value  of  the  minimum  CQM-tour,  LB4,  is  also  a  valid  lower 


a  '_a  "  *  ‘  ‘  m’  mf  •'  a*  _  *  -  •*_  ^  J  •**_  ■*" -  **  —  ' 


.S' 


bound.  The  dynamic  programming  recursive  equations  can  be  written  as 
follows.  Let 


f^(Q,U)  =*  length  of  the  minimum  t-path  from  node  1  to  i  which  picks  up  a 
total  load  of  Q  and  having  a  load  U  at  node  1. 


Then  for  i  «  2,...,M;  t  -  3,...,M;  Q  -  0 . QT  and  U  -  ),..., C 


fJ(Q,U)  *  Minimum 


fj  +  dji  if  i  c  S 

fj"1(Q,U+qi)  +  dJX  if  i  CD 


where 


fi(Q,U) 


if  Q  -  U  -  q1  and  i  e  S 


otherwise 


Then  LB4  can  be  computed  as  follows : 


LB4  »  Minimum  {^((^,0)  +  d^} 
i  eD 


The  CQM-tour  assures  that  any  minimun  M-tour  will  visit  a 
combination  of  supply  nodes  such  that  a  total  load  of  will  be  picked 
up  and  the  capacity  will  never  be  exceeded.  Further,  by  setting  U  *  0  as 
stage  M  the  M-tour  is  forced  to  visit  a  combination  of  demand  nodes  such 
that  a  total  load  of  will  be  delivered. 


-  > V  V'  . 
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The  M-tours  obtained  by  each  of  the  lower  bounding  procedures  will 
be  connected  tours.  However  these  tours  may  contain  subcycles  which  may 


% 


be  repeated  several  times. 

Considering  the  various  formulations,  1M  >  LB2  and  LB4  >  LB3,  since 
both  state  variables  of  LB2  and  LB3  are  imbedded  in  calculating  LB4. 
However  the  calculations  for  LB4  would  potentially  require  more 
computational  effort  compared  to  LB2  or  LB3. 

Each  of  the  previous  lower  bounds  can  further  be  strengthened  by 
appropriate  Lagrangian  ascent  procedures  (for  example,  see  Bazaraa  and 
Goode  [8]). 


gy 
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4.  A  BRANCH  AND  BOUND  ALGORITHM  FOR  THE  SINGLE  VEHICLE  PICK-UP  AND 
DELIVERY  PROBLEM 

The  lower  bounding  schemes  discussed  in  Section  3  were  implemented 
in  a  branch  and  bound  algorithm  for  the  single  vehicle  PDP.  Due  to  the 
lack  of  standard  test  problems  in  the  literature,  randomly  generated  test 
problems  were  used  for  testing  the  performance  of  the  algorithm.  The 
specifics  of  the  algorithm,  test  problems  and  the  computational  results 
are  given  below. 


A . 1  The  Dynamic  Programming  Network 

The  M-tour  problem  given  by  (3.1)  through  (3.6)  is  a  shortest  path 


problem  over  an  acyclic  M-partite  network.  In  this  network,  arcs  are 


defined  by  the  x^  variables,  where  each  arc  connects  two  different  nodes 


at  two  consecutive  stages.  If  there  had  been  no  precedence  and  capacity 
constraints  on  the  tours,  each  node  at  any  stage  will  be  linked  to  every 
other  node  at  the  next  stage.  By  making  use  of  the  special 
characteristics  of  the  problem  this  network  can  be  reduced,  thus  lower 
bounds  can  be  strengthened.  The  reductions  are  as  follows. 

For  convenience  suppose 


d^j  -  d^  for  i  “  j  -  t  «  1 ,  • . .  ,M 

a)  d!^  *  -  t  ■  2 , . . . ,M;  j  =  2,...,M 

dj!  "  C  j 


The  tour  must  s^art  and  end  at  the  base  node  1. 


b)  dij  -  "  3£D 

After  node  one  the  vehicle  can  only  go  to  a  supply  node. 

c)  db(r).(r)  '  *  *  E  S 

At  stage  M  the  vehicle  must  be  at  one  of  the  demand  nodes. 


d)  db( r)a( r)  “  *  r  "  1»***»N;  t  “  1 . M 

A  supply  node  a(r)  cannot  be  visited  immediately  after  its 

associated  demand  node  b(r). 


e) 


da(r)b(s) 


r  ^  s  *  1)..«)N 


At  stage  2  a  demand  node  of  an  order  cannot  be  visted  after  visting 
the  supply  node  of  another  order. 

f)  db(r)b( s)  “  *  r*  8  " 

At  stage  three  two  demand  nodes  cannot  be  visited  one  after 
another . 


g) 


M-2 

a(r)a(s) 


00 


r  *  s  -  1, . . .,N 


At  stage  M-2  two  supply  nodes  cannot  be  visited. 


h)  d 


M-l 

a(r)b(s) 


r  «  s  -  1, 


,N 


24 


At  stage  M-l  a  supply  node  must  be  connected  to  only  Its  associated 


demand  node. 


1)  Capacity  related  reductions: 

dC  .  CD 

a(r)a(s) 

d*.  ...  .  if  q  +  q  >  C  r^s  -  1 . N;  t 

b(r)b(s)  -  »  T  ^s  * 

dt  -  00 

a(r)b(s) 


1,...,M 


Two  nodes  cannot  be  connected  if  the  sum  of  their  loads  is  greater 
than  the  capacity  of  the  vehicle. 

In  enforcing  the  total  load  constraints,  i.e.,  equation  (3.7),  the 
variable  Q  at  each  stage  can  only  assume  certain  values.  For  example,  at 
stage  two,  Q  cannot  be  equal  to  the  total  load,  since  it  is  not  possible 
to  pick-up  all  the  load  at  the  first  pick-up  node  visited.  Likewise  at 
stage  M,  the  total  pick-ups  must  equal  to  Q^,.  For  the  vehicle  there  are 
two  extreme  node  visit  sequences  in  terms  of  pick-ups.  The  first  one  is 
to  visit  nodes  in  such  a  sequence  that  a  minimum  amount  goods  are  picked- 
up.  The  other  extreme  is  to  visit  nodes  in  an  order  such  that  a  maximum 
amount  of  goods  are  picked-up.  Suppose  the  tasks  are  ordered  with 
respect  to  the  nondecreasing  order  of  loads,  i.e.. 


qm  <  qi2]  < . <  q[N] 

Then  the  sequence  which  picks-up  the  minimum  amount  of  load  is  the 
one  that  visits  supply  and  demand  nodes  of  each  task  with  respect  to  the 
above  ordering  of  tasks,  Figure  1. 
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1 - > a [  1  ] >b(l  ] — >a[2  ] - >b[2] 


>a[N] - >b[N] - >1 


■^11]  ^[1]  +q[2j  “**[2]  +q[N]  “q[NJ 


Figure  1.  Representation  of  a  Vehicle  Tour  Which  Pick-up  Minimum  Load 

A  node  sequence  which  picks-up  the  maximum  amount  of  load  would  be 
the  one  which  visits  all  the  supply  nodes  first  and  then  visits  the 
remaining  demand  nodes.  In  this  case  the  supply  nodes  are  visited  with 
respect  to  the  noninceasing  order  of  task  loads.  Figure  2. 


1 - >a[N] - >a[N-l  ]— > . — >a[N] - >b(  i) - >b(  j)— . . . - >1 

"^[N]  +q[N-l]  +q[l] 


Figure  2.  Representation  of  a  Vehicle  Tour  Which  Pick-ups  Maximum  Load 

With  the  above  properties,  for  each  stage  the  minimisn  and  maximum 
number  of  possible  load  pick-ups  could  be  identified.  Let  be 

these  minimua  and  maximums  respectively.  Then 
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^  ■  iqi»i 


*2t+i  "  ht 
t-i 


nfi'iN-0-1] 


for  t  -  1, . .  .,N 


for  t  -  2 . N  +  1 


u 


t 


for  t  -  N  +  2, . . .  ,M 


At  each  stage  t  the  variable  Q  can  assume  values  <  Q  <  u  . 

The  2-node  subcycles  could  be  eliminated  without  increasing  the 
state  space,  by  the  procedure  given  Christofides  et  al  [16].  It  is 
sufficient  to  identify  the  best  and  the  next  best  t-tours  for  every 
possible  value  of  state  variables  at  each  stage. 

Using  the  above  described  state  space  restrictions  a  computer  code 
which  computes  minimum  CQM-tours  was  developed  by  Kirca  [23].  By  proper 
definition  of  the  state  variables,  namely  C  and  Q,  any  of  the  lower 
bounds  discussed  can  be  computed  by  this  algorithm. 

4.2  Branch  and  Bound  Algorithm 

A  branch  and  bound  algorithm  was  developed  for  the  single  vehicle 
PDP.  The  algorithm  is  capable  of  utilizing  different  lower  bounding 
schemes.  The  procedure  creates  partial  tours  and  computes  lower  bounds 
for  each  of  those.  If  the  lower  bound  for  any  of  these  tours  is  greater 
than  the  value  of  the  incumbent  solution  (upper  bound),  this  partial  tour 
is  discarded  from  further  consideration.  In  order  not  to  repeat  the 
already  evaluated  partial  tours  a  systematic  way  of  creating  and 
evaluating  partial  tours  is  required.  One  such  procedure  is  summarized 
below. 
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Step  0.  Initialization: 


Let  L  -  1,  Xl  -  D,  T  -  M  +  1,  G^T)  =■  1 
P,  -  {G.  (T)  },  Z*  -  P*  -  0 


Compute  a  lower  bound,  Z^,  for  the  partial  tour  P^.  If  Z^  >  Z  stop,  no 
feasible  solution  exists.  Otherwise  go  to  Step  1. 

Step  1 .  Candidate  tour  selection: 

*  * 

Pick  a  partial  tour  P^,  1  <  £■  <  L,  with  Z^  <  Z  .  If  none  exists  stop,  P 

* 

is  the  optimal  tour  with  value  Z  .  Otherwise  go  to  Step  2. 

Step  2 .  Branching : 

Let  T  »  M  +  2  -  2  P  „1 

for  each  i  e  X  repeat  steps  (2.1)  through  (2.3)  and  at  end  go  to  Step  3. 
Step  2.1.  L  *  L  +  1 


GL(t)  -  Gt(t) 
Gl(T)  -  i 

XL-Xi"  {1} 
if  ieD 


t  -  T  +  1 , . . . ,M  +  1 


X^  ■  X^U  {source  node  of  i } 

PL  "  ptD©L(T)> 

Step  2 .2 .  Compute  a  lower  bound ,  Z  ,  for  the  partial  tour  P  .  If  Z  > 

L  L  L 

* 

Z  let  Z^  ■  ®,  and  go  to  Step  2.2.  Otherwise  go  to  Step  2.3. 

_  *  * 

Step  2.3.  If  T  r  2  to  to  Step  2.2.  Otherwise  let  Z  ■  ZT  and  P  »  P. . 

L  L 

Go  to  Step  2.1. 

Step  3 .  Let  Z.  •  B  and  go  to  Step  1. 


By  selection  of  the  lower  bounding  scheme  for  Steps  0  and  2.2,  and  the 
partial  tour  selection  rule  for  Step  1,  the  branch  and  bound  algorithm 
can  be  specified.  The  algorithm  developed  here  utilizs  a  "depth  first" 
strategy.  The  partial  tour  which  has  the  least  lower  bound  among  the 
highest  level  active  partial  tours  is  selected  for  further  branching. 

This  procedure  minimizes  computer  storage  requirements.  At  any  instance, 
not  more  than  M  active  branches  exist  in  the  search  tree.  Furthermore, 
this  strategy  reaches  good  quality  upper  bounds  faster  so  that  the 
"fathoming"  efficiency  at  Step  2.2  is  increased. 

The  upper  bound  is  updated  at  Step  2.3  whenever  a  complete  tour  of  M 
nodes,  with  a  value  less  than  the  value  of  the  Incumbent  solution  is 
reached.  In  the  test  code,  updating  is  performed  whenever  the  lower 
bounding  procedure  yields  a  feasible  M-tour  with  a  value  less  than  the 
value  of  the  incumbent . 

The  algorithm  is  capable  of  computing  lower  bounds  with  any  of  the 
three  M-tour  lower  bounding  schemes  described  in  Section  3.  The  lower 
bounds  were  computed  by  the  algorithm  given  earlier  in  this  section;  Two 
different  versions  of  the  algorithm  were  developed.  The  first  version 
called  MULSOL,  solves  a  new  dynamic  program  over  a  new  updated  network 
for  each  of  the  partial  tours  evaluated.  In  the  second  version,  called 
ONESOL,  the  dynamic  program  is  solved  only  once  at  the  root  node,  i.e., 
at  the  initialization  step.  Then  the  dynamic  program  functional  values, 
f^(  )8»  are  stored.  Those  stored  values  are  later  used  to  compute  lower 
bounds  for  each  of  the  partial  tours  generated  in  Steps  2  through  4. 

Those  bounds  are  computed  as  follows. 
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The  advantage  of  ONESOL  over  MULSOL  is  that  it  requires  less 
computational  time  since  the  dynamic  program  is  solved  only  once. 

However  the  quality  of  the  lower  bounds  are  expected  to  be  poorer 
compared  to  MULSOL.  Further,  in  MULSOL,  storage  requirements  are  less 
than  ONESOL  since  the  dynamic  program  functional  values  need  not  be 
stored. 

4.3  Computational  Results 

The  branch  and  bound  procedures  were  tested  on  randomly  generated 
PDPs.  For  each  problem,  nodes  were  placed  randomly  on  a  two  dimensional 
space  and  the  distances  between  nodes  are  specified  as  integer  Euclidla. 
For  a  given  number  of  tasks,  N,  and  a  vehicle  capacity  C,  2N  +  1  nodes 
were  specified.  Node  1  is  fixed  as  the  base  node  and  supply  and  demand 

nodes  for  each  task  are  determined  as  (1  +  r,  N  +  r  +  1)  for  r  *  1 . N. 

The  load  for  each  task  was  drawn  from  a  uniform  distribution  between  1 
and  C.  The  details  of  problem  generation  and  problems  tested  can  be 
found  in  reference  [23]. 

A  second  set  of  problems  with  time  constraints  were  also  tested. 

Time  windows  on  the  pick-ups  and  deliveries  were  imposed  on  each  of  the 
test  problems  discussed  above.  The  time  windows  were  defined  in  terms  of 
stages  as  suggested  by  Armstrong.  A  vehicle  tour  has  2N  +  1  stages 
(stops).  By  specifying  “Not  Earlier  than  Delivery”,  NED,  and  "Not  Later 
than  Delivery",  NLD,  stages  for  each  task  time,  constraints  are 
implicitly  defined. 

All  algorithms  were  all  coded  in  FORTRAN  IV,  compiled  in  FTN 
(0PT»2),  and  executed  on  a  CDC  CYBER  730. 


The  initial  test  results  with  the  Branch  and  Bound  algorithms 
suggest  that  the  ONESOL  routine  is  faster  than  MULSOL.  MULSOL  was 
applied  only  with  the  CQM-tour  bounding  scheme. 

The  CQM-tour  routine  was  not  applied  in  ONESOL  because  of  the  excessive 
storage  requirements  for  larger  problems. 

The  algorithms  are  applied  to  the  5,  6,  7,  and  8  task  problems,  for 
each  task  size,  five  different  random  problems  were  tested.  In  Tables 
1  and  2  performance  of  the  procedures  are  tabulated.  A  120  cpu.  sec. 

Table  1.  Computational  Results  with  the  M-tour  Bounds  for  Single 
Vehicle  Problems  (N  ■  5  and  N  ■  6) 


PROBLEM 

N(M) 

« 

#  of 
Branches 

cpu. 

sec . 

#  of 
Branches 

cpu. 

sec . 

#  of 
Branches 

cpu . 

sec . 

51 

11 

13,003 

223 

.626 

207 

.731 

52 

26 

20.379 

277 

.645 

123 

.606 

5(11) 

53 

13 

15.705 

461 

1.040 

409 

1.101 

54 

11 

5.996 

39 

.286 

24 

.341 

55 

18 

16.847 

408 

.827 

135 

.562 

61 

72 

108.765 

30222 

43.882 

10929 

16.049 

62 

10  81 

120. 0002 

1767 

2.579 

1440 

2.709 

6(13) 

63 

8 

14.505 

519 

1.000 

156 

.999 

64 

34 

120. 0002 

2154 

3.306 

1549 

3.167 

65 

591 

120. 0002 

518 

1.400 

712 

2.136 

(1)  Optimal  reached  at  termination 

(2)  Terminated  at  this  time  limit 


Table  2.  Computational  Results  with  the  M-tour  Bounds 
for  Single  Vehicle  Problems  (N+7  and  N»8) 


ONESOL 

CM-TOUR  QM-TOUR 


PROBLEM 

N(M)  # 

if  of 
Branches 

cpu. 

sec . 

if  of 
Branches 

cpu. 

sec 

71 

1427 

2.588 

1487 

3.867 

72 

41524 

65.207 

58081 

92.504 

7(15)  73 

1923 

3.445 

1015 

3.078 

74 

9765 

15.494 

4928 

9.474 

75 

213 

1.026 

468 

2.638 

81 

1678 

4.056 

- 

- 

82 

31379 

54.201 

- 

- 

8(17)  83 

72865 

120. 0001 

- 

- 

84 

68814 

120. 0001 

- 

- 

85 

13094 

21.804 

- 

- 

(1)  Terminated  at  this  time  limit 

time  limit  was  set  for  all  problems,  and  whenever  this  limit  was  reached 
the  procedure  was  terminated  with  the  best  solution  found  so  far. 

Computational  results  with  ths  test  problems  suggest  that  the 
weakness  In  lower  bounds  with  ONESOL  is  highly  compensated  by  the 
computational  speed.  Furthermore  the  difference  between  the  CM-tour  and 
QM-tour  bounding  schemes  was  not  apparent.  Therefore  a  preference 
between  routine  is  more  favorable  for  the  CM-tour  procedure  as  compared 
to  the  (Jl-tour  procedure. 
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As  the  number  of  tasks  increases  the  solution  effort  also  increases 
for  all  procedures.  This  fact  suggests  that  the  size  of  problems 
solvable  by  these  procedures  is  limited  due  to  the  computational  effort 
required.  Also  the  storage  requirements  increase  as  the  problem  size 
Increases.  This  is  another  factor  limiting  the  successful  implementation 
of  the  exact  B&B  procedures  for  the  single  vehicle  PDPs  with  large 
numbers  of  tasks.  The  computational  effort  required  and  computer  storage 
space  requirements  for  different  schemes  are  tabluated  in  Table  3. 


Table  3.  Computational  Complexity  and  Memory  Requirements 
with  Different  Bounding  Schems  and  B&B  Strategies 


MULSOL  ONESOL 

D.P.  Branches  Memory  D.P.  Branches  Memory 

Routine  Effort  Effort  Space  Effort  Effort  Space 


M 

M2CQ 

C QM-tour 

M3CQ 

M3CQ 

2MCQ 

M3CQ 

QM-tour 

m3q 

M3Q 

2MQ 

m3q 

M 

m2q 

CM- tour 

m3c 

m3c 

2MC 

m3c 

m 

m2c 

Computational 

results  for  the 

time  constrained 

PDPs 

are  summarized 

Table  4. 

On  the 

unconstrained 

problems 

,  time  windows 

were 

Imposed . 
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Table  4. 

Time 

Constrained 
the  CM- tour 

Random 

Lower 

PDPs  Solved 
Bounds 

with 

PROBLEM 

N 

0 

Optimal 

Value 

Initial 

LB 

No.  of 
Branches 

Total  cpu. 
sec . 

51 

544 

427 

14 

.463 

52 

476 

472 

7 

.462 

5 

53 

520 

351 

90 

.535 

54 

629 

627 

8 

.392 

55 

346 

346 

0 

.398 

61 

529 

354 

125 

.852 

62 

551 

454 

29 

.664 

6 

63 

502 

19 

.638 

64 

664 

486 

52 

.823 

65 

428 

405 

44 

.721 

71 

518 

516 

7 

.819 

72 

603 

547 

207 

1.081 

7 

73 

699 

622 

46 

.908 

'4 


713 


672 


52 


85' 


The  effect  of  the  time  constraints  are  very  significant  in  the  sense  that 
they  reduce  the  state  space  for  the  lower  bound  computations,  thus  higher 
values  of  lower  bounds  were  reached.  Also  the  feasible  space  for  the 
overall  problems  is  reduced,  and  thus  the  search  effort  in  the  B&B  phase 
is  reduced.  As  a  result,  larger  size  PDPs  are  solvable  in  the  presence 
of  these  constraints. 

4.4  Bounds  Through  Lagrangian  Relaxation 

The  major  determinant  of  the  efficiency  of  the  B&B  procedure  is  the 
strength  of  the  lower  bounds  used.  As  the  size  of  the  problem  increases 
the  quality  of  the  lower  bounds  decreases  rapidly.  One  reason  the 
strength  of  the  lower  bounds  diminishs  rapidly  is  that,  as  the  problem 
size  increases,  the  subcycle  elimination  procedures  described  earlier 
loose  their  effectiveness.  Therefore,  the  M-tours  find  contain  subcycles 
which  repeat  themselves  several  times.  In  the  M-tours,  the  requirement 
to  visit  each  node  exactly  once  was  relaxed.  Hence,  the  procedure  does 
not  have  control  over  the  subcycles  or  force  the  tour  to  visit  each  node 
exactly  once.  One  way  to  enforce  these  node  visit  constraints  would  be 
to  dualize  them  and  apply  appropriate  Lagrangian  procedures. 

Suppose  we  relax  the  precedence  constraints  in  (2.8)  and  dualize  the 
node  visit  constraints  in  (2.5).  Then  the  resulting  program  is 

M  M  M  M  MM 

Max  F(w)  »  Min  LIE  d  x*  +  E  w  (  E  E  x*  -1  )  (4.1) 

t-1  i-1  j-1  1J  i-1  1  t-1  j-1  3 


s.t.  (2.2),  (2.3),  (2.4),  (2.6),  (2.7)  and  (2.9) 


The  objective  function  in  (4.1)  can  be  rewritten  as 

M  M  M  M 

Max  F(w)  ■  Min  £  £  £  (d.  ,  +  w  )  x*j ,  -  £  w  (4.2) 

w  t-1  i-1  j-1  J  J  i-1 

For  any  value  of  the  Lagrangian  multipliers  w,  the  solution  of  the 
above  program  in  (4.2)  results  in  a  CM-tour  problem,  where  the  tour 
visits  exactly  M  nodes  and  the  vehicle  capacity  constraints  are  not 
violated.  However  the  resulting  tour  may  visit  some  nodes  more  than  once 
and  some  nodes  may  not  be  visited  at  all. 

It  is  obvious  that  any  value  of  F(w)  is  a  lower  bound  for  the 
overall  PDP.  However,  larger  values  of  F(w)  will  yield  larger  lower 
bounds  which  could  help  significantly  to  reduce  the  search  effort  in  the 
B&B  phase.  For  this  purpose  a  subgradient  optimization  procedure, 
suggested  by  Bazaraa  and  Goode  for  the  TSP,  was  utilized  to  improve  the 
value  of  F(w).  The  algorithm  is  called  the  "Ascent  Procedure"  is  this 
discussion. 

The  algorithm  developed  for  the  single  vehicle  PDP  is  as  follows. 

STEP  1.  Apply  a  Lagrangian  Ascent  Procedure.  If  the  procedure 

★ 

results  in  a  feasible  PDP  tour  stop.  F  is  the  optimal  value  of  the  PDP 

* 

and  x  is  the  optimal  tour.  Otherwise  to  go  Step  2. 

STEP  2.  Branch  and  Bound  Procedure:  Apply  the  B&B  procedure 
described  in  Section  4.2  (ONESOL  version).  Replace, 
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for  all  l,j 


a.  .  «-  d.  .  +  w* 
ij  lj 


fJ(U)  -  f*C(U) 


for  all  i,t  and  0  <  U  <  C 


where  f^Uj's  are  the  dynamic  program  functional  values  correspondong 

*  * 

the  dual  and  primal  variables  (w  ,x  )  of  the  ascent  phase. 

The  procedure  was  tested  on  the  same  randomly  generated  problems 
given  in  Tables  5  and  6.  The  computational  results  are  encouraging. 


1  ^ 1 


-  ••  • . 
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Table  5.  Computational  experience  with  Single  Vehicle  Ascent+CM-tour 
routine  for  5  <  N  <  8 


(all  times 

are  in  cpu. 

sec . ' s) 

4  * 

■4  , 

PROBLEM 

Before 

After 

Optimal 

Ascent 

B&B 

Total 

k.  * 

N 

# 

Ascent 

Ascent 

Value 

Time 

Time 

Time 

LB 

LB 

\ 

51 

338 

446 

446 

2.291 

- 

2.291 

\  • 

52 

306 

452 

452 

1.254 

- 

1.254 

- 

5 

53 

337 

454 

454 

4.655 

- 

4.655 

:  : 

54 

510 

557 

557 

.474 

- 

.474 

55 

210 

305 

305 

2.709 

- 

2.709 

ie1 

61 

171 

451 

451 

7.742 

- 

7.742 

62 

356 

532 

533 

6.422 

.035 

6.457 

k  ,*  - 

6 

63 

319 

444 

444 

3.385 

- 

3.385 

i 

64 

453 

617 

617 

4.493 

- 

4.493 

■ 

65 

297 

317 

348 

6.506 

.073 

6.579 

‘‘‘.H 

.  *j 


■ 


Table  6.  Computational  Experience  with  Single  Vehicle  Ascent+CM-tour 


Routine 

for  9  <  N  < 

12 

(all  times 

are  in  cpu. 

sec .' s) 

PROBLEM 

N  # 

Before 

Ascent 

LB 

After 

Ascent 

LB 

Optimal 

Value 

Ascent 

Time 

B&B 

Time 

Total 

Time 

91 

269 

492 

537 

80.796 

5.365 

86.161 

92 

283 

356 

426 

74.557 

36.658 

111.215 

9  93 

525 

741 

748 

23.685 

.524 

24.209 

94 

394 

537 

597 

49.708 

17.704 

67.452 

95 

409 

589 

638 

55.554 

3.005 

58.559 

101 

334 

569 

620 

78.016 

5.312 

83.328 

102 

346 

576 

646 

89.089 

150.912 

240. 0001 

10  103 

457 

645 

680 

65.840 

2.864 

68.704 

104 

353 

498 

647 

102.705 

137.303 

240. 0001 

105 

306 

465 

546 

106.937 

19.488 

126.425 

121 

451 

592 

646 

232.983 

20.208 

253.191 

122 

365 

781 

826 

197.360 

40.312 

237.672 

12  123 

477 

673 

734 

140.545 

65.055 

205.600 

124 

449 

831 

868 

162.417 

26.455 

188.872 

125 

410 

655 

695 

105.568 

10.288 

115.856 

(1)  termination  at  this  time  limit 


are  encouraging.  In  Table  6  the  computational  results  with  the  above 
parameters  are  tabulated.  Column  (2)  gives  the  lower  bounds  obtained 
before  the  ascent  procedure.  Column  (3)  gives  the  lower  bounds  obtained 
at  the  termination  of  the  ascent  procedure. 
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As  can  be  seen  in  Table  5,  the  ascent  procedure  results  in  either 
the  optimal  solution  for  the  PDP,  or  produces  lower  bounds  so  close  to 
the  value  of  the  optimal  PDP,  that  the  B&B  phase  quickly  finds  and 
verifies  the  optimal  solution.  Encouraged  by  these  results  for  the  small 
problems  (5  <  N  <  8),  the  algorithm  was  tested  on  9,  10,  and  12  task 
problems.  The  results  are  given  in  Table  6,  for  randomly  generated 
problems . 

The  computational  results  suggest  that  the  Lagrangian  ascent  procedure 
produces  good  quality  lower  bounds,  at  the  expense  of  additional 
computational  effort.  Throughout  these  experiments  the  ascent  procedure 
resulted  in  good  quality  lower  bounds  (i.e.,  90Z  of  optimal)  in  early 
iterations.  Then  either  no  Improvement  was  achieved  or  the  Improvements 
were  not  significant  enough  to  justify  the  remaining  computational 
effort.  By  proper  selection  of  search  parameters  good  quality  lower 
bounds  can  be  obtained  more  efficiently.  Since  the  computational  effort 
in  the  B&B  phase  is  O(M)  at  each  branch,  a  trade-off  between  the  good 
quality  lower  bounds  and  B&B  procedure  could  be  investigated. 

Although  precedence  constraints  of  (2.8)  are  relaxed  for  the  above 
tests,  the  resulting  CM-tours  at  the  termination  of  the  ascent  phase  were 
either  feasible  or  near  feasible  with  respect  to  these  constraints.  This 
may  have  resulted  from  the  fact  that  in  order  for  the  M-tour  be  feasible 
with  respect  to  the  enforced  constraints  it  should  also  perform  pick-ups 
and  deliveries  such  that  the  precedence  constraints  are  satisfied,  (or 
least  violated) . 

In  each  of  the  ascent  iterations  a  CM-tour  is  solved.  The 

3 

computational  requirement  for  this  problem  is  0(m  C).  In  the  B&B  phase, 
the  lower  bounds  are  computed  by  utilizing  the  dynamic  programming 
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functional  values  found  at  the  end  of  the  ascent  phase.  In  this  case  the 
computational  complexity  at  each  branch  is  0(M).  The  computer  memory 
requirement  for  the  whole  procedure  is  0(^0  . 

4.6  Bounds  Through  1-trees 

If  the  capacity  and  precedence  constraints  of  the  PDP  model  of  Section 
2.1  are  relaxed,  the  resulting  problem  is  a  TSP.  Thus  any  optimal 
solution  or  a  lower  bound  for  the  TSP  is  also  a  lower  bound  for  the  PDP. 
Utilizing  TSP  lower  bounds  and  enforcing  the  capacity  and  precedence 
constraints  in  the  B&B  phase  could  be  one  approach  solving  the  PDP.  If, 
on  average,  the  value  of  the  optimal  PDP  tour  is  not  very  far  away  from 
the  value  of  the  optimal  TSP  tour,  then  1-tree  bounds  coupled  with  a 
duality  approach  may  result  in  good  quality  lower  bounds  for  the  PDP. 


4.7  Comparison  of  the  Exact  Single  Vehicle  Procedures 


¥ 

l 

Table  7.  Comparison  of  the  Computational  Effort  for  the 
Single  Vehicle  Methods 

(all  times  are  in  cpu.  sec.) 


Problem 

Size 

N 

CM- tour 
Total 
Time 

Ascent+( CM- tour) 
Ascent  Total 

Time  Time 

Ascent+(l-tree) 
Ascent  Total 

Time  Time 

5 

.689 

2.277 

2.277 

.167 

.548 

6 

10.430 

5.710 

5.717 

.240 

.671 

7 

17.552 

14.825 

14.825 

.392 

3.416 

8 

64.012 

38.138 

39.810 

.570 

9.293 

9 

- 

56.86 

69.609 

.571 

97.901 

10 

- 

88.517 

151.691 

- 

- 

12 

_ 

_ 

200.238 

— 

— 

Table  8.  Average  Strength  of  Lower  Bounds  Computed  by  Different  Schemes 

(lower  bound/ optimal  value) 


Problem 

Size 

CM- tour 

Ascent 

CM- tour 

Ascent 

1-Tree 

5 

.768 

1.000 

.917 

6 

.667 

.989 

.939 

7 

.665 

.987 

.881 

8 

.672 

.950 

.868 

9 

.638 

.921 

.807 

10 

.572 

.877 

- 

12 

.571 

.937 

— 
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clearly  that  the  procedures  which  utilize  Lagrangian  relaxation  are 
superior  to  the  one  that  does  not.  This  results  from  the  fact  that 
Lagrangian  relaxation  yields  good  quality  lower  bounds.  Comparison  of 
the  Lagrangian  CM-tour  and  1-tree  approaches  would  be  misleading.  By 
more  efficient  ascent  strategies  a  trade-off  between  the  asent  and  B&B 
phases  could  be  reached  for  the  CM-tour  scheme.  This  may  result  in 
reducing  the  computational  effort  at  the  ascent  phase  significantly.  As 
can  be  seen  for  the  eight  and  nine  task  problems  the  lower  bounds 
obtained  by  the  1-tree  approach  are  weaker  than  the  CM-tour  bounds.  Thus 
the  1-tree  routine  requires  more  computational  effort.  It  appears  that 
for  larger  problems  1-tree  bounds  will  be  weaker  and  thus  require  more 
computational  effort  compared  to  the  CM-tour  scheme. 

Thus  far,  in  the  literature  there  exist  only  two  exact  procedures 
for  the  single  vehicle  PDP.  These  are  the  procedures  due  to  Psaraftis 
[25]  and  Armstrong  [4].  Both  of  these  procedures  utilize  direct 
application  of  dynamic  programming  technique  on  the  PDP's.  In  Tables  9 
and  10  the  performance  of  these  two  methods  are  compared  with  the  ones 
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Table  9.  Comparison  of  the  Single  Vehicle  Procedures 
with  the  Methods  Exist  in  the  Literature 


Problem 

Size 

N(M) 


Existing  Methods 
Psaraftis  Armstrong 


Proposed  Methods(3) 

CM- tour  Ascent  Ascent 

CM-tour  1-tree 


!? 

2(5) 

.6 

3.92 

- 

- 

- 

3(7) 

.7 

4.08 

- 

- 

- 

2- 

4(9) 

1.1 

- 

- 

- 

- 

5(11) 

2.7 

5.23 

.685 

2.277 

.548 

"" 

6(13) 

9.6 

- 

10.430 

5.717 

.671 

■ 

7(15) 

46.8 

124.72 

17.552 

14.825 

3.416 

r 

8(17) 

149.7 

- 

64.012 

39.810 

9.293 

.  * 

9(19) 

591.4 

- 

- 

69.609 

97.901 

> 

(i) 

VAX /VMS  of 

MIT 

i 

(2) 

Includes  I/O  operations 

and  data  generation  on 

IBM  370/3303 

(3) 

Average  of 

five  random 

problems 

for  each  N  on 

CDC  Cyber  730 

Table  10.  Computer  Memory  Requirements  for  Different 
Single  Vehicle  PDP  Methods 


Problem 

Psaraftis 

CM-tour 

Ascent 

Ascent 

Size 

Armstrong 

CM-tour 

1-tree 

N 

2N3n-1+1 

(2N+1)2C 

(2N+1)2C 

( 2N+1 ) 

5 

811 

1210 

1210 

121 

10 

393611 

4410 

4410 

441 

15 

1.43  xlO8 

9610 

9610 

961 

10  20 

4.64  xlO 

16810 

1681 

30 

4.11xl015 

36000 

16810 

36000 


LC- 
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proposed  here.  Although  both  authors  report  their  computational  results 
for  the  dial-a-ride  problems,  for  the  same  number  of  tasks,  N,  the 
problems  they  report  and  the  problems  tested  here  may  be  considered 
essentially  equivalent. 

Due  to  differences  in  problem  assumptions,  computers,  and  coding,  a 
direct  comparison  of  the  three  methods  is  difficult.  Hence,  tables  9  and 
10  "suggest"  that,  for  the  instances  compared,  the  procedures  proposed 
here  seem  to  be  more  efficient  than  the  ones  in  Psaraftis  and  Armstrong. 
Even  the  least  efficient  method,  the  CM-tour  without  the  ascent  phase, 
performs  better  than  those  procedures  compared.  In  the  case  of  computer 
memory  requirements  the  procedures  suggested  here  (see  Table  8)  compare 
favorably  to  the  ones  available  in  the  literature.  In  the  algorithms  of 
Psaraftis  and  Armstrong,  memory  requirements  are  exponential  whereas  here 
the  proposed  methods  require  polynomial  memory  with  respect  to  the 
problem  size  N. 


5.0  CONCLUSIONS 
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